A method for direct observation of quantum tunneling in a single molecule. 
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f***-) , An application of impedance measurement technique (IMT) for a detection of quantum tunneling 

■ in molecular structures is investigated. A charged particle which tunnels in a two-well potential 

£\j ' is electrically coupled to a high-quality superconducting LC-circuit(tank) that makes possible a 

measurement of the electric susceptibility of the molecule at the resonant frequency of the tank. 
The real part of this susceptibility bears information about the tunneling rate through a measurable 
parameter - a phase angle between the tank voltage and a bias current applied to the tank. It is 
shown that the present approach is highly sensitive and allows the monitoring of the tunnel motion 
of charged nuclei in a single molecule. 
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Tunneling of a particle between minima of a double-well potential represents a pure quantum phenomenon which 
takes place in many physical and chemical systems 1 ' 2 . Splittings of molecular energy levels due to tunneling motion 
of constituent nuclei have been observed using spectroscopical methods 3,4,5 as well as methods of nuclear magnetic 
resonance relaxometry 6 . The measurements have been performed with samples containing a large number of molecules. 
Recent progress in molecular manipulation with a scanning tunneling microscope tip 7,8 makes possible a creation of 
single- molecule devices, and, in particular, apparatuses which use tunneling effects 9 . 

Here we propose to employ an impedance measurement technique 10,11 for contactless characterization of quantum 
tunneling in single molecules. This technique has provided a simple and reliable method for experimental investigation 
of tunneling in magnetic systems, especially, in superconducting flux qubits 12,13,14 . In the framework of IMT method 
the qubit's loop is inductively coupled to a high-quality LC circuit (tank) driven by a time-dependent bias current. 
An interaction of the tank with the qubit modifies an effective inductance of the tank and produces a shift of its 
resonant frequency. As a result, a voltage in the LC-circuit, (Vr(i)) — Vt cos(ujt + 0), has been displaced in phase 
from the bias current, Iuasif) — Iaccosut, by the angle 6. A dependence of the angle 9 on the bias applied to the 
qubit has demonstrated IMT dips which are indicative of quantum tunneling in the qubit. A tunneling rate between 
wells is determined by the width of the dips. From the theoretical point of view the tangent of the difference between 
phases of voltage and current, tan0, is proportional to the real part of a magnetic susceptibility of the qubit. 

To study intra-molecular tunneling of a charged particle, say, a proton or another light nucleus, in a two-well 
potential U(z), we have to measure an electric susceptibility of such a molecule. To do this, the molecule is placed 
between plates of the tank's capacitor Ct, so that the charged particle tunnels along lines of the intra-capacitor 
electric field E z — Vt / d. Here Vt is a voltage applied to capacitor's plates which are separated by the distance d. We 
> assume that the two-well potential U(z) has minima at z — ±zq. Taking into account the lowest tunneling doublet of 
' energy eigenstates we describe quantum dynamics of the system by the Hamiltonian 
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H = ja x + E -g z - (XV T + Q + f)o z + H T . (1) 
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Here A is a rate related to the particle's tunneling between the wells, e is a bias describing an asymmetry of the wells. 
This bias can be changed by applying an additional constant voltage to the capacitor Ct- A position of the particle 
is characterized by the Pauli matrix a z : z = z a z , so that an interaction with the intra-capacitor electric field E z 
is given by the term — XVt<t z , where A = e(zo/d). We add a qubit's internal heat bath with a variable Qo and an 
auxiliary external force /. It should be noted that thermal fluctuations of the tank voltage contribute to dccohcrence 
and relaxation of the qubit together with internal mechanisms. We consider the tank, having a resonant frequency 



lot = 1/\/LtCt, as a quantum harmonic oscillator with the Hamiltonian (h = 1): 



_<3 | H T = uj T (a + a + 1/2) - (a + a+)Q b - L T I T I bias + H TB . (2) 

The tank's excitations are described by creation/annihilation operators a + ,o([a, & + ]_ = 1); in so doing for the 
operators of voltage and current we obtain the expressions: Vt — iy / L0T/2Cr(a + — a), It = V W /2Lr(a + + a). A 
linewidth broadening jt and a finite quality factor of the tank, Qt = ujt/(^t), are caused by the tank's own heat 
bath, which is characterized by the variable Qb and the Hamiltonian Htb- 

Following to Ref.15 we derive the equation for the averaged voltage in the tank coupled to the double dot: 



2 



The particle's variable (a z ) is functionally-dependent on the tank's voltage: 



W *i<^y><M*i)>. (4) 



where the functional derivative corresponds to the electric susceptibility Xz( w ) °f the particle in the double-dot system 



The bias current, Iu a s{t) = -^ac c °swi, applied to the tank, induces oscillations of the tank's voltage, (Vr(t)) 
Vt cos(ujt + 9), with the amplitude Vr and the phase 6 which obey the equation 



-lo 2 - «7 T cj + u) T - '—^Xzi^) 
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V Te -ie = ^ Iac . (6) 

Or 



For the case when the frequency of the bias current is exactly equal to the initial resonant frequency of the tank, 
lo = lot, the current-voltage phase shift is given by the equation 

tane = -2A 2 ^ X ^ T ), (7) 
Or 

where Qt — wt/(2tt) is an effective quality factor of the tank, modified because of the double dot contribution to 
the broadening of the tank's linewidth: jt = 7t + {^ 2 oot / Ct)x"{oot)- 

The electric susceptibility of the double dot, Xz{oo), can De found in the framework of the theory of open quantum 
systems 15 ' 16 . Quantum dynamics of the charged particle in two-well potential is governed by the Heisenberg equations: 

a x = -sa y + 2(Q + f)<Ty, 

a v = -Aa z + e<r x - 2(Q + f)a x , 

a z = A<j y . (8) 

Here Q — Qo + XVr is the operator of the total dissipative environment which surrounds the double dot. For the 
susceptibility of this environment we have the relation: x( w ) — Xo( w ) + Xt(^), where X o(w) corresponds to the 
internal mechanisms of dissipation in the double dot, and 
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describes dissipative properties of the tank which are due to its coupling to the bath Qb- The spectral functions of the 
heat bath fluctuations, S(co) = Sq(u) + St(oj), is related to the imaginary part of the heat bath susceptibility by the 
fluctuation-dissipation theorem: S(lo) — x"{oo) coth(w/2T), with T being a temperature of the bath. Averaging the 
Heisenberg equations (8) over the equilibrium fluctuations of environment followed by an application of Bloch-Rcdfield 
approximation allows us to find the electric susceptibility of the charged particle confined in the two-well potential 

X.M = ~2A , , -™ + Si» + u*) + S(u,-uo) o (10) 

+ 7o)[-*(w - wo) + 7] [-t(w + lo ) + 7] 



Here lo = \/A 2 + e 2 is the ground-state splitting for the tunneling particle, a x = — (A/wo) tanh(wo/2T) is a steady- 
state value of the matrix a x . Coupling of the double dot to the dissipative environment results in the frequency- 
dependent decoherence and relaxation rates 7 and 70 

7 (w) - ^S(co) + 4 (l - — ) S(u> + wo) + 4 f 1 + — ) S(co - loo), 
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-f (co) = —[S(lo + lo ) + S(co -L0 )}. (11) 
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In particular, a back-action of the tank on the tunneling particle leads to the small decoherence and relaxation rates 

I / A \ ' .. e 2 /z n \ 2 / A • 



^,-i, ( 0,- l± l±) £ . (12) 
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It should be emphasized that the current- voltage angle 9 (7) is proportional to the double-dot susceptibility, ~x! z (u>t), 
taken at the resonant frequency of the tank lot, which is much lower than the tunneling frequency loq. Taking this 
fact into account we find the real part of the electric susceptibility: x'z^t) = — 2Act 9 /cj 2 . As a result, the phase angle 
between the voltage in the tank and the bias current is determined by the following equation 



By way of illustration, we consider tunneling in an ammonia molecule NH3 which is described by the two-well 
potential energy 17 . The existence of two delocalized states with eigenenergies separated by a tunneling splitting 
A/h — 23.8 GHz has motivated to harness this molecule as a quantum bit 9 . According to this proposal the NH 3 - 
molecule can be encapsulated into a fullerene Cqo to be addressed individually. A direct chemisorption of the ammonia 
molecule on the surface of a semiconductor represents another possiblity 18 . The fullerene ball incorporating the single 
ammonia molecule is inserted between the capacitor's plates having an area about 1 /jm 2 . If the distance between 
plates is of order of a diameter of Cqq— molecule 19 : d = 1.2 x 10~ 9 m , then for the capacitance Ct we get an estimation: 
Ct — 10 fF. Together with a value for the tank's inductance, Lt — 1 ^H, it gives the resonant frequency of the tank: 
wt/27t = 1.6 GHz which is much less than the tunneling frequency A/h : Ulut/A — 0.06. The electric dipole moment 
of the ammonia molecule 18 , 2ez , which is involved in the formula (13), is of order of 0.31 e x 10~ 10 m . For the case 
when losses in the tank are mainly due to a shunt resistance Rt, which is parallel to the capacitor Ct, the linewidth 
of the tank is defined as: jt = V(-RtCt), so that the quality factor is Qt = ljtRtCt/2, and the current-voltage 
angle 9 does not depend on the capacitance of the tank Ct ■ 



hu T Rt /2zqV A 3 , / \7~A~ 2 



tane = -^A {— ) (A 2+£2 )3/ 2 tanh ^— 2T~ J' (W) 

The current-voltage angle 9 can be measured in the experiment as a function of the bias e applied to the tunneling 
particle. With a reasonable estimation for the shunt resistance 19 : Rt = 50 x 10 6 Ohm, the single ammonia molecule 
demonstrates a pronounced low-frequency response: tanO = —0.23, at temperatures T < A/ks = 1.1 K, and e = 0. 
Measurement-induced decoherencc, 7 = 7(A) (12), is negligibly small in the process: 7/A ~ 10~ n , that allows us 
to characterize the IMT procedure as a weak continuous measurement 20,21 . To be realistic, for the example under 
discussion - the ammonia molecule confined inside of Cqq - we have to take into account screening effects due to a 
fullerene cage 9 . Because of the screening, a magnitude of the IMT dip (13,14) decreases that can be compensated by 
a proper increasing of the shunt resistance Rt and the quality factor of the tank Q T - 

In conclusion, we have proposed to use the impedance measurement technique for continuous monitoring of quantum 
tunneling in molecular systems which are characterized by a two-well potential energy. To this end, the molecule(s) 
should be inserted between the plates of a capacitor Ct comprising a high-quality superconducting tank along with 
an inductance Lt- We have shown that the present approach is sensitive enough to detect a tunneling motion in a 
single molecule. 
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